in 
o 
o 

(N 

> 
O 

in 



> 

00 
O 



in 
o 



X 



Mesons at high temperature in Nf = 2 QCD* 
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We report first results for spectral functions of charmonium in 2-flavour QCD. The spectral functions are 
determined from vector and pseudoscalar correlators on a dynamical, anisotropic lattice. J/ip and 77c are found to 
survive well into the deconfined phase before melting away at T < 2Tc. Current systematic uncertainties prevent 
us from drawing any definite conclusions at this stage. 



1. INTRODUCTION 

The properties of hadrons or hadronic reso- 
nances above the deconfinement transition is a 
subject at the heart of the current experimen- 
tal programme at RHIC, where hadronic signals 
are used to obtain information about the state of 
matter inside the fireball. The questions of in- 
terest include the issue of which hadrons survive 
as bound states in the quark-gluon plasma, and 
up to which temperature; as well as the trans- 
port properties of light and heavy quarks in the 
plasma. 

These properties are all encoded in the spec- 
tral functions p{uj,p), which are related to the 
imaginary-time correlator Gr (t, p) according to 



1 

Gt{t,p)^— pr{uj,p)K{T,ui)dui 
Jo 



(1) 



where the subscript F correspond to the different 
quantum numbers. The kernel K is given by 



KiT,u) 



cosh[cj(r - 1/2T)] 

sinh(w/2T) (2) 
e""nsH + e-""[l + nBH]. 



and nB is the Bose-Einstein distribution func- 
tion. 

The spectral function can be extracted from 
the lattice correlators G{t) using the Maximum 
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Entropy Method (MEM) [T]. For this to work 
and give reliable results, it is necessary to have a 
sufficient number of points in the euclidean time 
direction. This will be prohibitively expensive, 
especially when dynamical quarks are included, 
unless anisotropic lattices are used, with a tempo- 
ral lattice spacing much smaller than the spatial 
lattice spacing. 

We will here focus on the charmonium S-wave 
states rjc and J/^ at zero momentum, which have 
attracted much attention following the suggestion 
121 that J/ip suppression could be a probe of de- 
confinement. Potential model calculations using 
the heavy quark free energy have tended to sup- 
port this picture. However, previous simulations 
in the quenched approximation |8l4l5j indicate 
that contrary to this, J/ip may survive up to tem- 
peratures as high as 1.5 — 2Tc. Recently, potential 
model calculations using the internal energy of 
the heavy-quark pair have reached the same con- 
clusion, and using the most recent lattice data 
|S] these models indicate a qualitatively similar 
picture in the case of Nf = 2 QCD |6I7| . 

In this study we attempt to determine directly 
the spectral functions of charmonium in 2-flavour 
QCD using anisotropic lattices and the Maximum 
Entropy Method. 

2. SIMULATION DETAILS 

We use the Two-plaquette Symanzik Improved 
gauge action 8 and the fine- Wilson, coarse- 
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Hamber-Wu fermion action W with stout-link 
smearing JO]. The process of tuning the ac- 
tion parameters, and the parameters used, are de- 
scribed in more detail in |11I12) . The parameters 
correspond to a spatial lattice spacing as ~ 0.2fm 
with an anisotropy ^ — as/ at ~ 6. The sea quark 
mass corresponds to m,r/wp ~ 0.55. 

These actions are designed for large 
anisotropics, and a quenched study found 
that no mass-dependent tuning of the quark 
anisotropy was necessary up to valence quark 
masses well beyond charm. This appears no 
longer to be the case for dynamical quarks; in- 
deed, for the parameters used here the anisotropy 
determined from the pion dispersion relation was 
found to be ^ 6.4, in rough agreement with 
the gluon anisotropy, while the anisotropy from 
the charmonium dispersion relation was found to 
be ~ 8. This issue is still under investigation. 

For this preliminary study we have used a spa- 
tial volume of = 8'^ and generated 100 config- 
urations, sampled every 10 HMC trajectories, for 
Nt = 48, 32, 24 and 16. Nt = 32 corresponds to 
a temperature close to the pseudocritical temper- 
ature Tc, although larger lattices will be needed 
to determine this with any precision. We have 
computed charmonium correlators in the pseu- 
doscalar (jjc) and vector {J/tp) channels with bare 
charm quark mass atuic = 0.1. The charmonium 
spectrum in the hadronic phase is presented in 
|12I13| . In this study we have used local (un- 
smeared) operators, 

Gr(r) = '^^4i^^ t)Mr{y, t + t)) , (3) 

x,y,t 

where 

Mr(f , t) = i,{x, t)r^j{x, t), r - 75, J^ . (4) 

All-to-all propagators ^1] have been used to im- 
prove the signal and sample information from the 
entire lattice. The propagators were constructed 
with no eigenvectors and two noise vectors diluted 
in time, colour and even/odd in space. 

The MEM analysis has been performed with 
the continuum free spectral function as default 
model, using the euclidean correlators in a time 
window starting at r = 2, and cutting off the 
energy integral in ^ at atLOmax = 6. 
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Figure 1 . Average Polyakov loop as a function of 
temperature. 



3. RESULTS 

In an attempt to locate the pseudocritical tem- 
perature on these lattices, we performed simula- 
tions varying Nt in the range 28-40, measuring 
the real part of the Polyakov loop (L) and its 
susceptibility. The results for the Polyakov loop 
are shown in fig. ^ As can be seen from the 
figure, the Polyakov line in the transition region 
follows a linear behaviour in temperature, and 
therefore no pseudocritical temperature can be 
determined on these lattices. It will be necessary 
to use larger lattices to determine Tc- For this 
reason, and because of the uncertainties regard- 
ing the anisotropy, we will refrain from quoting 
results in terms of Tc- 

The free-fermion spectral functions are shown 
in fig. 121 for Nt — 32 and 24, and in the con- 
tinuum. The most striking feature of the lattice 
free spectral functions is the spike at atw ~ 0.6. 
This means that lattice artefacts are big at this 
point, and results for spectral functions cannot 
be trusted in this region. The lattice functions 
undershooting the continuum curve appears to 
be a mass-dependent, 0{atm) or ©(a^m), ef- 
fect, evidenced by the effect being much larger 
at atm = 0.2 [T^. 

The spectral functions obtained from the MEM 
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Figure 2. Free-fermion spectral functions in the 
pseudoscalar and vector channels, for bare quark 
mass at Too = 0.1. 



Figure 4. Vector {J/tp) spectral function for dif- 
ferent temperatures. 
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Figure 3. Pseudoscalar (r\c) spectral function for 
different temperatures. 



the deconfinement temperature, finally melting 
away at T < IT,^. This is in qualitative agree- 
ment with recent potential model results using 
the static quark-antiquark internal energy [H]. 
Using the colour singlet internal energy \J\ yields 
a dissociation temperature of Tdis ^ 2Tc, while 
a potential constructed to exclude the gluon in- 
ternal energy [TB] yields Tdis ~ ^AT^ 0. Given 
the uncertainties in this calculation, both these 
results are consistent with the present data. 

It is not clear whether the apparent stronger 
binding in r]^. for iVt = 24 is significant, in par- 
ticular given that no such effect is observed for 
J/t/i. Higher statistics will be needed to resolve 
this. If it is confirmed, it might be in line with 
potential model calculations using J7i, which is 
strongly peaked around Tc- 



analysis are shown in figs. |31 and ^ for ric and J/i/; 
respectively. The main peak position for Nt = 48 
agrees within errors with the mass obtained for 
the respective particles on the same lattices, using 
a variational basis of smeared operators l[I2J,I^ . 
The second peak coincides with the cusp in the 
free spectral function in fig. 13 and it can there- 
fore be concluded that this is primarily a lattice 
artefact. The radially excited states cannot be 
resolved with our present statistics. 

The results in figs.|niand0]indicate that the IS 
states survive in the medium up to well beyond 



4. OUTLOOK 

Using an anisotropic lattice, we have performed 
the first calculation of charmonium spectral func- 
tions in 2-flavour QCD. The results appear to 
confirm the picture emerging from quenched sim- 
ulations, that J/'0 and r]c survive until well into 
the deconfined phase. 

These simulations have been performed with 
parameters that are not fully tuned, giving rise to 
significant systematic uncertainties. A particular 
issue is the tuning of the charm quark anisotropy. 
We are in the process of obtaining a fully tuned 
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parameter set, and simulations at these parame- 
ters will be carried out in the very near future. 
Firm conclusions will have to await these simula- 
tions. 

The small lattice volume — only (1.6fm)'^ — 
is also a major source of systematic uncertainty, 
in particular as it has prevented a determination 
of the pseudocritical temperature. Future simula- 
tions will be carried out on larger spatial volumes. 
It will also be important to increase the statistics, 
in order to resolve excited states and disentangle 
the effects of thermal width and finite statistics. 
We will also be carrying out a systematic study 
of the effects of using different default models, in- 
cluding the free lattice spectral functions, and dif- 
ferent time and energy ranges in the MEM anal- 
ysis. Initial indications are that our results are 
relatively robust against such changes. 

The coarse spatial lattice is an issue in that it 
gives rise to significant, mass-dependent, lattice 
artefacts as shown in fig. |2] It will ultimately be 
necessary to repeat the calculation on a finer lat- 
tice; this will however require a new nonperturba- 
tive tuning and is therefore not on the immediate 
horizon. Some information about lattice spac- 
ing effects may be gleaned from quenched stud- 
ies using the same action, which are considerably 
cheaper as the quark and gluon anisotropics can 
be tuned independently. 

We are planning to compute the spectral func- 
tions also at non-zero momentum 117! , which con- 
tain additional information not found at zero mo- 
mentum, and which may relate more directly to 
experimental data taken at non-zero momentum. 

We will also study light vector meson corre- 
lators at zero and non-zero momentum, which 
can be related to dilepton production rates in the 
plasma. Finite lattice spacing effects are expected 
to be less of a problem in this case, so the current 
lattices are likely to be suitable for this purpose. 
This will proceed once the fully tuned parameter 
set is ready. 
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